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Abstract 

We present a lattice calculation of the B meson i?-parameter Bb using the 
NRQCD action. The heavy quark mass dependence is explicitly studied over 
a mass range between mi, and 4mb with the ©(l/mg) and ©(l/mg) ac- 
tions. We find that the ratios of lattice matrix elements (Ojy*)/(^o'*)^ 
(O^*) / (Aq"*)^, which contribute to Bb through mixing, have significant I/tuq 
dependence while that of the leading operator (0^*)/(Aq"*)^ has little l/mq 
effect. The combined result for BB{mb) has small but non-zero mass depen- 
dence, and the BB{mi,) becomes smaller by 10% with the l/mg correction 
compared to the static result. Our result in the quenched approximation at 
/?=5.9 is Bb^{5 GeV) = 0.75(3)(12), where the first error is statistical and 
the second is a systematic uncertainty. 

PACS number(s): 12.38.Gc, 12.39.Hg, 13.20.He, 14.40.Nd 
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I. INTRODUCTION 



The constraints on the unitarity triangle of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix can provide us with one of the most crucial information on the physics beyond the 
Standard Model [|I|]. However, due to large theoretical or experimental uncertainties, the 
current bound is too loose to test the Standard Model or new physics. The B^-B^ mixing 
sets a constraint on iV^t^Vj^l through the currently available experimental data on the mass 
difference between two neutral B mesons AM^= 0.477±0.017 ps~^ 0. The experimental 
achievement is impressive as the error is already quite small ~ 4%. Theoretical calculation 
to relate AM^ to |VJdV^fe|, on the other hand, involves a large uncertainty in the B meson 
matrix element {B'^\0 i\B^) , which requires a method to calculate the non-perturbative QCD 
effects. 

Lattice QCD is an ideal tool to compute such non-perturbative quantities from first prin- 
ciples. There has been a number of calculations of the B meson decay constant fs and the 
S-parameter Bb {Bb describes the matrix element through {B^\Ol\B^) = (8/3)i?s/^M^). 
The calculation of /b is already matured at least in the quenched approximation Ma- 
jor systematic errors are removed by introducing the non-relativistic effective actions and 
by improving the action and currents. Remaining a (lattice spacing) dependent systematic 
error is confirmed to be small, and in some papers the continuum extrapolation is made. 
A recent review summarized the value of J'b as /s= 165±20 MeV within the quenched 
approximation. 

An essential ingredient of these calculations is the use of the non-relativistic effective 
actions. Since the 6-quark mass in lattice unit arrib is large for the typical lattices used for 
simulations, the relativistic (Wilson-type) actions could suffer from a large discretization 
error of order arrib or {amb}'^. The non-relativistic QCD (NRQCD) [Q, on the other hand, 
is formulated as an expansion in p/rriQ. In the heavy-light meson system, where the typical 
spatial momentum scale is A.qcd, the error from the truncation of higher order terms is 
controllable. The calculation of fs is now available to order 1/ml, and it is known that the 
contribution of 0{l/mb) is significant (~ —20%) while that of 0(l/mf) is small (~ —3%) [Q. 
In addition, the calculations based on the Fermilab approach for heavy quark |^,[^ agree 
with the NRQCD results including their l/rriQ dependence. These results make us 
confident about the non-relativistic effective action approaches in Lattice QCD. Now that 
the computation of Jb is established, the next goal is to apply the similar technique to the 
computation of Bb- 

The lattice calculations of the 5-parameter have been done in the infinitely heavy quark 
mass (static) limit [P-|Tl[], and the results are in reasonable agreement with each other. 
There is, however, some indication that the l/rrib correction would be non- negligible from 
the study with relativistic quark actions |T^Jl3[] . Their results show that there is small but 
non- zero negative slope in the l/mg dependence of 5-parameter, but it is not conclusive 
because of the possible systematic uncertainty associated with the relativistic quark action 
for heavy quark. The purpose of this work is to study the l/rriQ dependence of Bb by 
explicitly calculating it with the NRQCD action at several values of l/mq. Our result 
confirms the previous works |ll2| , [l3[| : there is a small negative slope in Bb- In addition, 
we find that the slope comes from the large l/mq dependence of and -B^**, which are 
matrix elements of non-leading operators. For the observed l/mq dependence of the lattice 
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matrix elements, qualitative explanations are given in Discussion section using the vacuum 
saturation approximation. 

The perturbative matching of the continuum and lattice operators introduces a compli- 
cation to the analysis. Since the one-loop coefficients for four-quark operators are not known 
yet for the NRQCD action, we use the coefficients in the static limit in Refs. |T^-|18|. The 



systematic error associated with this approximation is of a^j [amo) and expected to be small 
compared to the l/mQ correction itself. This and other systematic errors are discussed in 
detail in the Discussion section. 

This paper is organized as follows. In the next section, we summarize our matching 
procedure. The simulation method is described in section |T|, and our analysis and results are 
presented in section The results are compared with the vacuum saturation approximation 
on the lattice in section [V|, and we estimate the remaining uncertainties in section [V^. Finally 
our conclusion is given in section [VI I|. An early result of this work was presented in Ref. [jl9| . 



II. PERTURBATIVE MATCHING 



In this section, we give our notations and describe the perturbative matching procedure. 
The mass difference between two neutral Bq mesons is given by 



-6/23 



where q = d or s and S(){xt) {xt = rn^/m^) and are so-called Inami-Lim function and the 
short distance QCD correction, respectively. Their explicit forms can be found in Ref. 
is a Ai?=2 operator 



(2) 



renormalized in the MS scheme with the naive dimensional regularization (NDR). is 
related to the anomalous dimension at the next-to-leading order with nj active flavors as 



Jrif 
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where 



/3o 



11 - 3^/, 
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A = 102 - ^Uf, 
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-7+ -Hf. 
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(4) 



nf=5 when /i is greater than or equal to the b quark mass. 
The S-parameter Bb„ is defined through 



(5) 



where denotes the axial-vector current b^^'^bq- The renormalization invariant B- 
parameter is defined by 
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Bf 
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(6) 



which does not depend on the arbitrary scale fi up to the next-to-leading order. The scale 
/i is conventionally set at the scale of 6-quark mass fi = nib. 

In order to calculate the matrix element {B^\OL{mb)\B^) on the lattice, we have to 
connect the operator (^^(mf,) defined in the continuum renormalization scheme with its 
lattice counterpart. The matching coefficients can be obtained by requiring the perturbative 
quark scattering amplitudes at certain momentum with continuum Ol operator and with 
lattice four fermi operators should give identical results. 

At the one-loop level the matching gives the following relation 



OLirrib) 



a. 



l + _±[41n(aX: 



Dl - 14] ) O'til/a 



where C^{i*ij7vs} denotes the naive local operators defined on the lattice in which the light 
quarks are not rotated. Their explicit forms are the following, 

Or = &7m(1 + 75)g&7M(l + 75)g, 

On = &7m(1 - 75)'?&7m(1 + 75)g + ^7^(1 + 75)g&7M(l - 75)'? 
+26(1 - 75)g6(l + 75)g + 26(1 + 75)^6(1 - 75)?, 

Os = 6(1 - 75)g&(l - 75)g. (7) 

Unfortunately, the one-loop coefficients -D{L,i?,Ar} for the NRQCD heavy and 0(a)-improved 
light quark action pl[ are not known yet. In this work, we use the one-loop coefficients in 



the static limit pSHTgl, 

Dl = -21.16, 



D 



R 



-0.52, 



D 



N 



-6.16, 



(8) 



The systematic error associated with this approximation is at most as/ {airiQ), since the 
NRQCD action's m ^ 00 limit agrees with the static action. The numerical size of the error 
is discussed later. 

The matching of the axial-vector current appearing in the denominator of Eq. (^ can 



be done in a similar manner 15,22,23 



= (1 + ^[21n(aX) + Da- \]) <*(l/c 



(9) 



where Aq"* is defined on the lattice, and the matching coefficient Da in the static limit 
Da = -13.89. 

In calculating the ratio of Eq. (^ a large cancellation of perturbative matching corrections 
takes place between the numerator and denominator, since the large wave function renormal- 
ization coming from the tadpole contribution in the lattice theory is the same. To make this 
cancellation explicit we consider the matching of a ratio Bsirrib) = (Cl("^6))/(8/3)(Ao)^ 
itself. 
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+Z^/A2K; l/a)5j^*(l/a) + Z^/^^K; l/a)5^''*(l/a), (10) 

where 5{f ,R,^,5}(l/a) = (Og,i?,w,5}(l/«))/(8/3)(^[r*(l/a))^ and the 5 and 5 states are 
understood for the expectation values (■ ■ ■) as in Eq. (|^). Then the coefficients become 

Z^/A.K; 1/a) = (^1 + ^(D^ - 2D A - f )) , (11) 

OL 

ZR/A^inib, 1/a) = -^Dr, (12) 
ZN/A^imb, 1/a) = -^Dn, (13) 
Z5/A2(mfc;l/a) = ^D5. (14) 
The Eqs. (pJ] ) -([l^ ) are used in the following analysis to obtain BB{mb)- Numerical values 

of 

Z{L,R,N,s}/A^{''T^b] 1/a) are given in Table | for the lattice parameters in our simulation. 
For the coupling constant as in Eqs. (|TTD-(p!4D we use the V-scheme coupling p4[ with 
q* = n/a or q* = 1/a. At /3=5.9 those are av/(7r/a)=0.164 and ay(l/a)=0.270. The tadpole 
improvement does not make any effect on the ratio of Eq. ([l^), since the tadpole 



contribution cancels between the numerator and denominator. The 6-quark mass scale rrib 
is set to 5 GeV as usual. 

In the previous works in the static approximation P-[TT[], the leading and the next-to- 



leading logarithmic corrections are resummed to achieve better control in the running from 
rrih to 1/a. In this paper we use the one-loop formula without the resummation for simplicity. 
This does not introduce significant error, since the mass scale difference between rrib and 1/a 
is small and the effect of the resummation is not important. In Appendix ^ we compare the 
Z factors with and without the resummation. 

We determine the heavy-light pseudo-scalar meson mass Mp from the binding energy 
^bin mgagu];gcl in the simulation using a formula 

aMp = Z^aniQ - Eq + aE^'"", (15) 

where Z^ and Eq are the renormalization constant for the kinetic mass and the energy 



shift, respectively. Both have been calculated perturbatively by Davies and Thacker ||25 
and by Morningstar [^. Since the precise form of their NRQCD action is slightly different 
from ours, we performed the perturbative calculations for our action. Our results for the 
coefficient A and B in the perturbative expansion 



Eo = aviq*)A, (16) 
Zrn = l + av{q*)B, (17) 



are summarized in Table |I|. 



III. SIMULATION DETAILS 



Our task is to compute the ratios {O'-^l j^j^gy{l/ a)) / (8/3) {A^Q^{l/a))'^ using lattice 
NRQCD with the lattice spacing a. In this section we describe our simulation method 
to obtain them. 
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We performed the numerical simulation on a 16^ x 48 lattice at (3= 5.9, for which the 
inverse lattice spacing fixed with the string tension is l/a= 1.64 GeV. In the quenched 
approximation we use 250 gauge configurations, each separated by 2,000 pseudo-heat bath 
sweeps. For the light quark we use the 0(a)-improved action |^ at /t=0.1350, 0.1365. The 



clover coefficient is set to be Csw = l/^o, where uq = (Ppiag)^/'^=0.8734 dX (3 = 5.9. The 
critical k value is =0.1401, and corresponding to the strange quark mass determined 
from the K meson mass is =0.1385. 

For the heavy quark and anti-quark we use the lattice NRQCD action with the tadpole 
improvement — > U^/uq. The precise form of the action is the same as the one we used 
in the previous work We use both of the ©(l/mg) and 0{l/m^) actions in parallel in 
order to see the effect of the higher order contributions. The heavy (anti-)quark field in the 
relativistic four- component spinor form is constructed with the inverse Foldy-Wouthuysen- 
Tani (FWT) transformation defined at the tadpole improved tree level as in Ref. f^. The 
heavy quark masses and the stabilization parameters are (amg, n) = (10. 0,2), (5.0,2), (3.0,2), 
(2.6,2) and (2.1,3). These parameters approximately cover a mass scale between 4mb and 
rrib. 

We label the time axis of our lattice as t = [—24, 23] . The heavy quark and anti-quark 
propagators are created from a local source located at the origin (t=0 on our lattice) and 
evolve into opposite temporal directions. The light quark propagator is also solved with 
the same source location and with a Dirichlet boundary condition at t = —24 and t = 23. 
The B and B mesons are constructed with local sink operators. Thus we have the four- 
quark operators at the origin and extract the matrix elements from the following three-point 
correlation function 



c^5f^(^i,'^2) = EE(oK*^'^i,^i)o!,"*(o,oX*^(t2,f2)|o), (is) 

XI X2 



where X denotes L, R, N or S. Because of a symmetry under parity transformation, 
C^^\ti,t2) and C^^\ti,t2) should exactly coincide in infinitely large statistics. Therefore we 
explicitly average them before the fitting procedure we describe below. 

To obtain the ratios i?^*(l/a) we also define the following two-point functions 

c^'\ti) = E(oK*(ti,:?X*^o,o)|o), (19) 

X 

c^'\h) = E(oK*(o,oX*^(t2,f)|o), (20) 

X 

and consider a ratio 



|C(2)(ti)C(2)(t2) |(P0|^[,»*(l/a)|0)(0|4'^*(l/a)|P0) 



Bf{l/a), (21) 



where P° denotes a heavy-light pseudo-scalar meson. The ground state P° meson is achieved 
in the large \ti\ regime. Although we use the local operator for the sinks at ti and t2, the 
ground state extraction is rather easier for finite amg than in the static approximation, 
since the statistical error is much smaller for NRQCD p7| , p8| . This is another advantage of 
introducing the l/rriQ correction. 
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The physical BBijnh) is obtained by extrapolating and interpolating each i?^*(l/a) to 
the physical B meson with k and mg, respectively before combining them as Eq. ([To|). The 
final result for i?s(mft) may also be obtained by combining the ratio of correlation functions 
before a constant fit. Namely we use the relation 

Bp{mb-MM) = Zx/A^'^bA/o) ^ f \r{l)(+ \ ~^ Bp{nn,). (22) 

X=L,R,N,S S*-^ '\^2) 

Since the statistical fiuctuation in the individual _B^*(l/a) is correlated, the error is expected 
to be smaller with this method (We use the jackknife method for error estimation). Following 
Ref. |]ll| we refer to this method as the "combine-then-fit" method, while the usual one as 
Eq. ( pll) is called the "fit-then-combine" method in the rest of the paper. 



IV. SIMULATION RESULTS 



We describe the simulation results in this section. The results are from the 0{l/mQ) 
action unless we specifically mention. 



A. Heavy-light meson mass 



The binding energy of the heavy-light meson is obtained from a simultaneous fit of two 
two-point correlation functions. The numerical results are listed in Table |Tl| for each amq 
and K. Extrapolation of the light quark mass to the strange quark mass or to the chiral 
limit is performed assuming a linear dependence in 1/k. 

The meson mass is calculated using the perturbative expression Eq. ([151) . The results 
with Q;y(7r/a) and with av{^la) are also given in Table fT|. 



B. B^f{l/a) and Bp{mb) 



Figures |1| and ^ show the ti dependence of Bp{mb]ti,t2) in the "combine-then-fit" 
method. The perturbative matching of the continuum and lattice theory is done with the 
l^-scheme coupling ay(g*) at q*=7i/a (left) and 1/a (right). Their difference represents the 
effect of 0(0;^). The signal is rather noisier at amQ=5.0 (Fig. |l]) than at amQ=2.Q (Fig. ^ 
from the same reason as in the static limit P7|,|2B|. But, still, a reasonably good signal is 



observed even for large amq. A plateau in the ti dependence is reached around ti= 8 ~ 
11 for both ^2 = ~10 and —15. To be conservative we take as well as \t2\ greater than 
10 for the fitting region. All data points (^1,^2) in 10 < < 13 and in 10 < |t2| < 13 
are fitted with constant to obtain the result for Bp{mb)- We confirm that except for the 
heaviest quark the results are stable within about one standard deviation under a change 
of the fitting region by at most two ti steps in the forward and backward direction. The 
numerical results are listed in Table [TV] . 

The light quark mass (1/k) dependence of Bp{mb) is presented in Fig. |^. Since its 
dependence is quite modest, we assume a linear dependence on 1/k and extrapolate the 
results to the strange quark mass and to the chiral limit as shown in the plot. Results of 
the extrapolation are also listed in Table 1^. 
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C. 1/Mp dependence 



In Fig. ^ we plot Bp{mb) in the chiral limit, namely Bp^{mi,), as a function of 1/Mp in 
the physical unit. We take q*=n/a (circles) and 1/a (squares) for the scale of ay- Regardless 
of the choice of the coupling, we observe a small but non-zero negative slope in 1/Mp, which 



supports the previous results by Bernard, Blum and Soni [|T2[ and by Lellouch and Lin |T3 
using the relativistic fermions. 

To investigate the origin of the observed 1/Mp dependence, we look into the contributions 
of the individual operators 0''^l rn s} through the "fit-then-combine" method with the same 
fitting region as before. We list the results for each in Table |V[ Figure |^ shows the 
1/Mp dependence of S["*(l/a)(=Sj^"*(l/a)), B^^\l/a) and 5^"*(l/a). While no significant 
1/Mp dependence is observed in i?^*(l/a), i?j^*(l/a) and i?g*(l/a) have strong slope. Since 
their sign is opposite and the sign of the matching factors Z^/A^i^b] ^/a) and Zs/A'^if^b] 1/a) 
(see Table |) is the same, a partial cancellation takes place giving a small negative slope for 

We also make a comparison of the results of the l/mq action (circles) with those of the 
I/toq action (triangles) in Fig. |^. There is a small difference between the two results in 
B^§\l/a) and in 5(f*(l/a) toward large 1/Mp (1/Mp > 0.2 GeV"^), which is consistent 
with our expectation that the difference is an 0{h.QCD/^QY effect. 

Previous results in the static approximation by Ewing et al. (diamond) , Gimenez and 
Martinelli (triangle) [|10| and Christensen, Draper and McNeile (circle) |]ll| are plotted with 
filled symbols at 1/Mp = 0. Although the (3 value and the light quark action employed (the 
0(a)-improved action is used in Refs. [^|10| and unimproved action in Ref. [|ll[]) are different, 
all the results are in good agreement with each other. A quadratic extrapolation (dashed 
line) using our l/mg NRQCD result also does agree with these previous static results. 



D. Result for Bpimb) 

Combining the data for S^*(l/a) discussed above, we obtain Bp^{mf,) with the "fit-then- 
combine" method. We confirm that the difference in numerical results from both methods 
are completely negligible. Figure |^ shows the results of "fit-then-combine" method using 
av{l/a) with both actions. The comparison with the static results is also made in this plot, 
where only the statistical error in each calculation is considered and the same matching 
procedure as ours are applied. We again observe a consistent result. 

Interpolating the above NRQCD results to the physical B meson mass, we obtain the 
physical BpS'^b) 

R /^A_/0.r8(3) (g* = 7r/a) . . 

5B,(m,) - I Q_^2(3) (g* = l/a) ^^^> 

for the Oil/mo) action, and 

^^^("'"^ - I 0.71(3) (g* = l/a) (2^) 
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for the 0{l/mQ) action. The quoted error is statistical only. For the ratio of BbJBb^, we 
obtain BbJBb^ = 1.01(1) for q* = n/a and BbJBb^ = 1.02(1) for q* = 1/a from both 
actions. 



V. DISCUSSION 

The strong 1/Mp dependence in i?^*(l/a) observed in Fig. ^can be roughly understood 
using the vacuum saturation approximation (VSA) on the lattice as explained below. Here 
it should be noted that the terminology of VSA we use here does not immediately mean 
BB{mb) = 1. 

The VSA for is unity by construction. This is true even for finite 1/Mp, and its 
prediction is shown by a straight line in Fig. ^(a). The NRQCD data is located slightly 
below the line (~0.9), but the mass dependence is well reproduced by the VSA. 

For B'§^ and we require a little algebra to explain their mass dependence under the 
VSA. Using the Fierz transformation and inserting the vacuum, we obtain 

(po\o^at^pO^ = -^(i^|67M75g|0)(0|67M75g|P°) - y(i^|&75g|0)(0|675g|P°), (25) 
(po\0'at^p0^ = l{nh5q\0){0\h5q\P'), (26) 

where \P^) denotes a heavy-light pseudo-scalar meson at rest, and (0|67^75g|P'') is related 
to the pseudo-scalar decay constant 

{nA',^\l/a)\0) = (0K*(l/a)|p0) = fpMp. (27) 

Let us now consider a decomposition of the 6-quark field b into the two-component non- 
relativistic quark Q"^ and anti-quark x fields. Up to ©(l/mg) we have 

i^,, = (Qt o)(l + I-^),.,-(0 X){l + l^],.,, (28) 



and then 



{P'\b^,q\0) = (PO|(gT 0) I 1 + ^ I 75^10) 

(P^|&7075g|0) + 2(P0|(gt o)^75g|0), (30) 

(0|675g|P°) = -(o|(o x) 1 1 + 75g|P°) 

■(0|67o75g|P°) + 2(01(0 x)l^l,q\P')- (31) 
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By defining Sfp as 



2m, 



Q 



-(o|(o 



x)|^75g|P°) 



2m, 



SfpMp 



Q 



we obtain 



(po|Oj^*|P°) 
(po|0!."*|p0) 



^/|M^ (^1 - 
-IfpMl (l 




(32) 

(33) 
(34) 



In our previous work we denoted 5fp as (5/p 



(2) 



Thus the VSA for Pj^* 



and for P^* read 



B 



Zat{VSA) 



'7> 



B 



/at(VSA) 



(35) 
(36) 



neglecting the higher order contribution of order l/mg. Results for Sfp/fp is available at 
(3=5.8 in Ref. §. We plot them in Fig. |(b) and |(c) by crosses, which show a qualitative 
agreement with the measured values. 

Di Pierro and Sachrajda [l^ pointed out that the value of several P-parameter-like 
matrix elements of the B meson is explained by the VSA surprisingly well in the static limit. 
Here we find that the l/mg dependence of _B^*(l/a) can also be reproduced qualitatively. 
This result suggests that the vacuum saturation is a reasonable qualitative picture for the 
heavy-light meson. It does not mean, however, that the VSA works quantitatively for Bb{ij), 
and careful lattice studies are necessary for precise calculation of the S-parameters. 



VI. REMAINING UNCERTAINTIES AND THE FINAL RESULT 

To estimate the systematic uncertainties in lattice calculations is a difficult task. In our 
case it is even more true, since we have a simulation result only at a single (3 value. However 
we attempt to do it, giving a dimension counting of missing contributions. 

The following sources of systematic errors are possible: 

• discretization error: Both of the heavy and light quark actions are 0(a)-improved 
at tree level, and there is no discretization error of 0{aKQCD)- The leading error 
is of Oi^o^Is^Qfjj-)) and of 0{asaKQCD)- The second one is from the missing one- loop 
perturbative correction in the 0(a)-improvement (its matching coefficient has been 
already obtained in Ref. |T^). We naively estimate the size of them to be 0{a'^K^Qj^) ~ 
O^aKqcDCts) ~ 5%, assuming 1/a ~1.6 GeV, Aqcd ~300 MeV and ~0.3. 

• perturbative error: The operator matching of the continuum and lattice AB=2 oper- 
ators are done at one-loop level. Thus the 0(a;^) correction is another source of error. 
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In addition, we use the one-loop coefficient for the static lattice action, though our 
simulation has been done with the NRQCD action. The error in this mismatch is as 
large as 0{as/ {amo)) and 0{asIs^QCD/i^Q)- The size of these contributions is 0{al) 
~ 0{aJ{amQ)) ~ 10% and ©(a^AQc^Mg) ~ 2%. 

• relativistic error: Since we have performed a set of simulations with the ©(l/mg) 
action and the 0{l/rnQ) action, we can estimate the error in the truncation of the 
non-relativistic expansion. As we have shown, the difference between the results with 
0{l/mQ) and 0{l/mQ) is small (~ 2%) around the B meson mass. Then the higher 
order (©(l/mg)) effect is negligible. 

• chiral extrapolation: We have only two light quark k values. Then the linear be- 
havior in the chiral extrapolation is nothing but an assumption. Although the light 

quark mass dependence is small and the assumption is a reasonable one, we conserva- 
tively estimate the error from the difference between the data at our lightest k value 
(«;=0.1365) and Kc- It leads 3% for BB^irrib). 

• quenching error: All results are obtained in the quenched approximation. Study of 
the sea quark effect is left for future work. 

Taking them into account, we obtain the following values as our final results from the 
quenched lattice, 

Sb,K) = 0.75(3)(12), 
1^ = 1-01(1)(3), 

where the first error is statistical and the second a sum of all systematic errors in quadrature. 
In estimating the error of the ratio BbJBb^ we consider the error from chiral extrapolation 
only, assuming that other uncertainties cancel in the ratio. The above result is related to 
the scale invariant S-parameter Bb^ as 

= 1.12(4)(18) 



~ ^ Mm)]-'/'' 



BB^mb) = 1.15(5)(18) ' ^^^^ 



using the leading and next-to-leading formula, respectively, where we use AQf;,^=0.237 GeV 
and the two-loop /3-function. 



VII. CONCLUSION 

In this paper we investigate the (^(Aqcd/^^q) and ©(Ag^^/mQ) effects on the B- 
parameter. We find that there is no significant mass dependence in the leading operator 
contribution i?^"*(l/a), while the mixing contributions -B^*(l/a) and i?^"*(l/a) have large 
O^Kqcd/^q) corrections. The O^I^^Qj^/mq) correction for each 5^*(l/a) is, however, rea- 
sonably small for the B meson as we naively expected. The observed l/mq dependence is 
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qualitatively understood using the vacuum saturation approximation for the lattice matrix 
elements. 

The lattice NRQCD calculation predicts the small but non-zero negative slope in the 
mass dependence of Bp{mb) and about 10% reduction from static limit to the physical 
B meson. In the present analysis, we combine lattice simulation for finite heavy quark 
mass with the mass independent matching coefficients determined in the static limit. The 
dominant uncertainty is, therefore, arising from the finite mass effects in the pcrturbative 
matching coefficients. For more complete understanding of the I/ttiq dependence, matching 
coefficients with the finite heavy quark mass are necessary. 
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APPENDIX A: 



In this appendix, we compare our perturbative matching by simple one-loop formula with 
the renormalization group (RG) improved ones used in the previous static calculations 



Since the matching procedure for determining the RG improved coefficients is given in 
Refs. in detail (see also Refs. [p^ -|3T|), we just show the results appropriate for our 



actions and definition of operator. 

Considering the matching of a ratio BB{rnh) = {OL^rrih)) / (8/3) {Aq}"^ again as in sec- 
tion H, the RG improved versions of Zx/A'^{''T^b', l/o) are as follows, 



ZL/A^im,; 1/a) = Z^' ( 1 + ^(I^l - Wj,) ) , (Al) 



ol 

Zr/a^ {m,- 1/a) = Zr' X -^Dr, (A2) 



QJ 

Zn/a^ (m,; 1/a) = Z^' x -^D^, (A3) 
Zs,AAm,-l/a) = Zr\ (A4) 



and 



47r ' 3 
Oisijni,) If as{mb 



8/25 




+ ^(-8) I I 7^ ) - 1 1 7, (A5) 

= ^4^(-8) ( 1 . (A6) 

Numerical values of Z^ l,r,n,s} ia^ i^b] 1/*^) given in Table | together with those by the sim- 



ple one-loop formula, where we use the V-scheme coupling as appearing in Eqs. ([Al| 



while the couplings in Eqs. ( [A5|) and ( [Aq ) are defined in the continuum MS scheme 
with A^=344 MeV, which corresponds to A^= 237 MeV. 

MS ' ^ MS 

Now assuming each is of 0(1), the dominant effects of resummation arise from 
Zl/A"^ and Zs/A"^. Since, however, its difference is at most 5% level and the effects from 
Zi/A"^ and ZsjA? are destructive, the total effect amount to less than 3%. To be specific, 
using our data extrapolated to the static hmit (see Table |IV|) to calculate B^°-^{mh), we 
obtain the results tabulated in Table |V| from the two matching procedures. In this case the 
effect of resummation is almost negligible. 
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TABLES 



Q 


av{q*) 










One-loop 


ir/a 


0.164 


0.973 


-0.080 


-0.007 


-0.104 




1/a 


0.270 


0.956 


-0.132 


-0.011 


-0.172 


RG improved 


ir/a 


0.164 


0.930 


-0.069 


-0.006 


-0.118 




1/a 


0.270 


0.978 


-0.113 


-0.010 


-0.118 



TABLE I. Matching factors at /?=5.9 by the two different procedures. 



amq 


10.0 


5.0 


3.0 


2.6 


2.1 


A 


1.011 


0.946 


0.855 


0.819 


0.754 


B 


-0.075 


0.018 


0.119 


0.152 


0.329 



TABLE II. One-loop coefficients for the self-energy. 



amq 


10.0 


5.0 


3.0 


2.6 


2.1 


Ac=0.1350 
K=0.1365 


0.691(10) 
0.655(12) 


0.675(4) 
0.636(5) 


0.664(3) 
0.625(3) 


0.660(3) 
0.620(3) 


0.653(2) 
0.613(3) 


K,. =0.1 385 
Kc=0.1401 


0.608(15) 
0.571(18) 


0.586(7) 
0.547(8) 


0.574(4) 
0.534(5) 


0.569(4) 
0.529(4) 


0.560(3) 
0.520(4) 


Mp[GeV]{q*=TT/a) 
Mp[GeV]{q*=l/a) 


16.864(30) 
16.558(30) 


8.867(13) 
8.719(13) 


5.662(8) 
5.576(8) 


5.018(7) 
4.944(7) 


4.279(7) 
4.268(7) 



TABLE III. The numbers in second and third lines are the binding energies obtained from 
the simultaneous fits. Those in fourth and fifth lines are the results of linear extrapolations to 
the strange and chiral limit. The numbers in last two rows are the corresponding physical meson 
masses with g*=7r/a and 1/a. 
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amq 



B 



lot ^ 
lat I 



lat( 



(=^r(«)) 



Bpinib) 
with ay(7r/a) 



Bpinib) 
with Q!y(l/a) 



10.0 
5.0 

3.0 
2.6 
2.1 



0.94(5 
0.91(2 
0.92(1 
0.92(1 
0.92(1 



1.27(8) 
1.47(4) 

1.85(3) 
2.00(3) 
2.27(4) 



-0.69(3 
-0.73(1 

-0.82(1 
-0.86(1 
-0.93(1 



0.88(4 
0.84(2 
0.83(1 
0.82(1 
0.80(1 



0.84(5 
0.79(2 
0.77(1 

0.75(1 
0.72(1 



10.0 

5.0 

3.0 

2.6 
2.1 



0.95(6 
0.90(3 
0.91(2 
0.91(2 
0.91(1 



1.25(12) 
1.45( 5) 
1.85( 4) 
2.01( 4) 
2.30( 5) 



-0.69(4 
-0.73(1 
-0.82(1 
-0.86(1 
-0.93(1 



0.89(6 
0.83(3 
0.82(2 
0.81(2 
0.79(1 



0.85(6 
0.78(3 
0.76(2 
0.74(2 
0.71(2 



oo 
10.0 
5.0 
3.0 
2.6 
2.1 



0.93(11 
0.97( 9 
0.89( 4 

0.90( 2 
0.90( 2 
0.90( 2 



1.15(23) 
1.22(17) 
1.43( 7) 
1.86( 6) 
2.03( 5) 
2.33( 6) 



-0.68(7 
-0.69(5 
-0.72(2 

-0.82(2 
-0.86(2 
-0.93(2 



0.92(14 
0.91( 9 
0.82( 4 

0.80( 2 
0.79( 2 
0.78( 2 



0.89(14 
0.87( 9 
0.77( 4 
0.74( 2 
0.73( 2 
0.70( 2 



oo 
10.0 
5.0 
3.0 
2.6 
2.1 



0.94(14 
0.98(11 

0.88( 4 
0.89( 3 
0.89( 2 
0.89( 2 



1.11(28) 
1.20(21) 
1.42( 9) 
1.87( 7) 
2.04( 6) 
2.35( 7) 



-0.68(8 
-0.69(6 

-0.72(2 
-0.82(2 
-0.85(2 
-0.93(2 



0.94(17 
0.93(11 

0.81( 4 
0.79( 3 
0.78( 2 
0.77( 2 



0.91(17 
0.89(11 

0.76( 5 
0.73( 3 
0.72( 3 
0.69( 2 



TABLE IV. Numerical results for -Bf *(a) and Bp{mb) 





One-loop 
q* = ir/a q* = 1/a 


RG improved 
q* = tt/o q* = 1/a 




0.90 0.88 


0.88 0.88 



TABLE V. Numerical results for S^*"*(m(,). Statistical errors are omitted here. 
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FIG. 1. Bp{mi,]ti,t2) as a function of ti, while t2 is fixed at —10 (circles) or —15 (squares). 
The heavy quark mass is amQ=5.0 and /t=0.1365. The "combine-then-fit" method is used. In the 
perturbative matching, av{TT/a) and ay(l/a) are used in the left and right plot respectively. 
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FIG. 2. Same as Figure |l| but at amQ=2.6. 
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FIG. 3. Extrapolation of Bp{mb) to the strange and to the chiral hmit. The heavy quark mass 
is 0777,(5=5.0 (circles) and 2.6 (squares). Results with g*=7r/a (left) and 1/a (right) are shown. 
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FIG. 4. Inverse heavy-light meson mass dependence of Bpj{mb) with q*=-iT/a (circles) and 1/a 
(squares). 
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FIG. 5. Inverse heavy- light meson mass dependence of (a) -B^*^, (b) -Sjy^*, (c) B^g^. For what 
symbols and line denote, see text. 
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FIG. 6. Comparison of our NRQCD data (open symbols) with the static ones (filled symbols). 
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